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Abstract 

Domain wall structure which may form in theories with spontaneously broken parity is generically 
in conflict with standard cosmology. It has been argued that Planck scale suppressed effects can 
be sufficient for removing such domain walls. We study this possibility for three specific evolution 
scenarios for the domain walls, with evolution during radiation dominated era, during matter 
dominated era, and that accompanied by weak inflation. We determine the operators permitted by 
the supergravity formalism and find that the field content introduced to achieve desired spontaneous 
parity breaking makes possible Planck scale suppressed terms which can potentially remove the 
domain walls safely. However, the parity breaking scale, equivalently the majorana mass scale Mr 
of the right handed neutrino, does get constrained in some of the cases, notably for the matter 
dominated evolution case which would be generic to string theory inspired models giving rise to 
moduli fields. One left-right symmetric model with only triplets and bidoublets is found to be 
more constrainted than another admitting a gauge singlet. 
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I. INTRODUCTION 



Existence of right handed neutrino states [1, 2, 3, 4] is a strong indicator of parity as a 
fundamental symmetry of nature, spontaneously broken in low energy physics. The scale of 
parity breakdown is as yet unknown. The see-saw mechanism [5, 6, 7, 8] while providing 
elegant qualitative explanation is unable to make a narrow prediction of the relevant energy 
scale due to wide variation in the fermion masses across the generations. In this paper we 
study the role of domain walls formed by spontaneous parity breaking in determining or 
constraining the energy scale of breakdown of this symmetry. 

To be specific we study the implementation of parity in Left-Right symmetric models [9, 
10, 11, 12, 13]. While such models can descend from an SO(10) grand unified theory 
[14, 15, 16, 17, 18], the scale of such unification is known to be high(~ 10 16 GeV). On 
the other hand the scale of parity breaking is completely undetermined and could be much 
lower. Indeed there are no observational obstructions to the parity breaking scale, and 
the associated right handed neutrino majorana mass scale, from being as low as TeV scale 
[19, 20, 21]. Here we study a particular source of constraint on this scale imposed by 
cosmology, with the possibility of restricting the parity breaking scale to low values. It 
suffices for this purpose to focus on the Left-Right symmetric model alone, independent of 
how the model may unify into SO (10). For the purpose of protecting the low scale theory 
from large radiative corrections we impose supersymmetry. Specifically we investigate the 
constraints placed on models incorporating TeV scale supersymmetry. 

A robust consequence of approximate left-right symmetry in the early universe is the 
occurance of transient domain walls. It has been proposed [22] that such domain walls are 
susceptible to instability arising from non-renormalizable operators suppressed by Planck 
scale. Supergravity then introduces two interesting ingredients not considered in the previous 
treatments of domain walls. Firstly the structure of the non-renormalizable terms is dictated 
by the supersymmetry formalism[23, 24, 25]. On the other hand one has to contend with 
the danger of gravitino overabundance [26, 27]. In this paper we explore the restrictions on 
the possible energy scale of parity breaking imposed by these considerations. 

We study these effects in the context of two implementations of left-right symmetry, one 
where all superfields carry non-trivial gauge couplings and another, for comparison, which 
contains a gauge singlet. We also study these models within three different scenarios for 
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the dynamics of the wall complex. One is a scenario in which the walls disappear within 
the radiation dominated era, another where dominance of moduli keeps the universe matter 
dominated during the domain wall evolution, and a third wherein the domain walls in fact 
come to dominate the universe for a limited epoch, accompanied by a mild inflationary 
phase. In all the scenarios of domain wall evolution, the left-right symmetric model with 
a singlet turns out to be less restricted than the one without singlets. The model without 
any singlet turns out, at least in one scenario of wall evolution, to be sufficiently restrictive 
that the parity breaking scale can be no larger than 10 8 GeV. The overall lesson is that the 
new features introduced by supergravity can have a strong bearing on the scale of parity 
breaking for ensuring viable cosmology free of permanent domain walls. 

The alternative to supergravity induced terms for distinguishing between parity symmet- 
ric vacua was studied in [28, 29] wherein the parity breaking operators are induced at a much 
lower scale, viz., the supersymmetry breaking scale, and signalled by the gauge mediation 
mechanism, thus linking the two scales. This avenue for evading unwanted domain walls 
remains open for models that get restricted in the present study. 

In the remainder of the paper, in sec. II we review cosmology with domain walls, in sec. 
Ill we discuss three possible scenarios for evolution of domain walls and the requirement 
in each case for the successful destabilization of the wall complex. In sec. IV we discuss 
the origin of the parity breaking terms within supergravity formalism. In sec. V we discuss 
the essentials of supersymmetric left-right model. In subsection V A we discuss a particular 
renormalisable implementation of left-right supersymmetric model and then check for the 
sufficiency of the induced Planck scale terms to cause the required destabilisation of the 
domain walls. The same study is carried out for a recent implementation of supersymmetric 
left-right symmetry including a gauge singlet in subsection VB. In sec. VI we summarise 
the conclusions. Two appendices, appendix A and appendix B contain the calculations of 
the Planck scale suppressed terms for the two left-right symmetric models studied. 

II. DISCRETE SYMMETRY BREAKING AND COSMOLOGY 

Spontaneous breakdown of discrete symmetries in unified models can give rise to domain 
walls. Stable domain walls from unified theories have long been recognized as inconsistent 
with the observed Universe [30]. In the presence of several degenerate ground states the 
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domain wall network in the Universe can be rather complicated. If the domain walls are 
stable, the energy density stored in the network decreases as p oc 1/a, resulting in a(t) oc t 2 
leading to a mild inflationary behavior. A more generic possibility is that the walls con- 
tinue to be destroyed due to collisions and result in the formation of homogeneous domains. 
However even one domain wall of grand unified scale of the size of the horizon can conflict 
with cosmology. For domain walls arising at symmetry breaking scale Mr, it can be esti- 
mated [30, 31] that the density perturbation introduced by them would conflict with known 
magnitude of temperature fluctuation AT/T ~ 10~ 5 of the cosmic microwave background 
if Mr > lMeV. This impasse is overcome if the spontaneously broken discrete symmetry is 
also broken explicitly by a small amount. For example, [32] the symmetry <p — » —<\> can be 
broken by adding a term e0 3 to the Lagrangian which gives a pressure difference governed 
by the small parameter e between the two sides of the domain walls. 

The authors of [22] have discussed several similar reasons for considering such gravity 
induced terms and their effect in destabilizing domain walls. For the theory of a generic 
neutral scalar field <j), the effective higher dimensional operators can be written as 

Such terms give rise to a pressure difference across a given domain wall of the amount of 
the difference in the effective energy density across the wall, 5p = AV e //. Specifically, the 
terms odd in <j) break the discrete symmetry, and the leading contribution to the difference 
in pressure is fj^(0) 5 - From cosmo logical considerations we can separately estimate the 
difference 5p in the energy density across a domain wall required for timely removal of the 
domain walls. It is found that this has a value smaller by a factor Mr/Mpi than the leading 
order term in the generic effective potential considered above, leading to the conclusion that 
the walls will indeed be removed without conflicting with cosmology. 

This toy example however is only instructional because in realistic theories, gauge invari- 
ance and supersymmetry significantly constrain the structure of terms that can arise. In a 
non-supersymmetric example [33], gauge invariance implies that the leading order operator 
is suppressed by l/M Pl . Further, the terms are products of vacuum expectation values of 
different scalar fields which may differ significantly in their mass scales, as will be true in 
our study. The exceptional case where the toy example may be of direct relevance is the 
presence of one or more gauge singlet scalar fields in the theory, ( equivalently, superfields 



in a supersymmetric theory) permitting the kind of terms listed above. 

As an illustration of this phenomenon consider the leading order operator containing 
several scalar fields, 0j 

0102030405 / x 

Borrowing from a calculation that will be detailed later (see eq. (13) ), suppose the required 
constraint for successful removal domain walls is 5p > Mp 1 ^ 2 /M P { 2 . Then the requirement 
that the operator in Eq.(2) is sufficient for removing domain walls is that 



^102030405 > M R 



11/2 



M Pl ~ Mp{ 2 



(3) 



Now suppose that there are only two kinds of fields, one getting the vacuum expectation 
value of the order of Mr, the parity breaking scale, and that there are x factors of this field 
in the operator, while the other field constitutes the remaining factors, and gets a TeV scale 
value v. Then 

M 2 

M 1 
1vi pi 

so that, taking TeV scale to be v ~ 10 3 GeV, 

, „, / 24.5 - 3x 

\ogM R <— (5) 

5.5 — x 

This relation means that if x — 5, Mr can be as large as the Planck scale, while for x — 1, 
Mr is forced to have value < 10 5 GeV, 

In our analysis we shall be assuming that parity breaking occurs at the same mass scale 
as the mass of the heavy majorana neutrinos. The constraints to be derived also depend 
on a few ancillary details, specifically the dynamics of the walls before they disintegrate. 
The primary implication of these other details, to be discussed in the following, is only 
the value of the temprature at which the standard cosmology resumes. In the following 
we shall admit the possibility of this more general kind of evolution and focus attention 
on two issues. The Universe should be radiation dominated at temperature lOMeV and 
lower in order to ensure successful Big Bang Nucleosynthesis(BBN). Secondly, the danger 
of gravitino overabundance is generic to all supersymmetric models. Detailed calculations 
[26, 34] show that the gravitinos with unacceptable consequences to observable cosmology are 
generated entirely after reheating of the universe subsequent to primordial inflation provided 
Tr > 10 9 GeV. Thus we make the requirement that entropy generated from the decay of 
domain walls should not raise the temprature of the Universe above this temperature scale. 



III. MODELS OF DYNAMICS OF DOMAIN WALL COMPLEX 



Occurance of domain walls per se at some epoch in the early Universe is not ruled out, 
provided the walls eventually get destroyed. Safe disappearance of domain walls was dealt 
with in some generality in [22] and [35], the former in the context of Planck scale effects, and 
the latter in the context of instanton induced effects from QCD. There have been several 
model specific studies of the fate of domain walls, e.g., [36, 37, 38] and studies pointing 
out that transient domain walls may in fact form the basis for explanation of some of 
the cosmological effects such as leptogenesis, [33, 38, 39, 40] or address problems such as 
proliferation of relics[28, 37, 38]. 

For the purpose of this paper we consider three possible routes through which domain 
walls may evolve. The first one consists of domain walls originating in radiation dominated 
universe and destabilized and destroyed also within the radiation dominated era before 
they begin to dominate the energy density of the Universe, referred to as RD model. The 
second scenario was essentially proposed in [37], which consists of the walls originating in a 
radiation dominated phase, subsequent to which the universe enters a "matter dominated" 
phase, either due to substantial production of heavy unwanted relics such as moduli, or 
due to the presence of a coherently oscillating scalar field. This we refer to as MD model. 
Here too walls disappear before they come to dominate the energy density of the Universe. 
Finally, we consider a variant of the MD model in which the domain walls come to dominate 
the energy density of the Universe and continue to do so for a considerable epoch, leading 
to mild inflationary behaviour or weak inflation [41], [42]. We refer to it as WI model. We 
now describe these in detail. 

A. Evolution in radiation dominated universe 

The essentials of this scenario are as originally proposed by Kibble [30] and Vilenkin [32] . 
Domain walls arise at some temperature T c , the critical temperature of a phase transition 
at which a scalar field <fi acquires a non-zero vacuum expectation value at a scale Mr. The 
energy density trapped per unit area of the wall is a ~ M\. The dynamics of the walls is 
determined by two quantities, force due to tension /j> and force due to friction fp. The first 
of these is determined by intrinsic energy per unit area <r, and the average scale of radius 
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of curvature R prevailing in the wall complex. We estimate fa ~ cr/R. The frictional force 
is proportional to the collisions encountered by the wall with surrounding radiation with 
energy density ~ T 4 , while the former is navigating through the medium at speed j3. This 
force is estimated as fp ~ j3T A . The epoch at which these two forces balance each other sets 
the time scale ~ R/(3. We may take this as the time scale by which the wall portions that 
started with radius of curvature scale R straighten out. Putting together these statements 
we get the following scaling law for the growth of the scale R(t) on which the wall complex 
is smoothed out. 

R(t) « (Ga) 1/2 t 3/2 (6) 

Now the energy density of the domain walls goes as pw ~ (aR 2 /R 3 ) ~ (a/Gt 3 ) 1 ^ 2 . In 
radiation dominated era this pw becomes comparable to the energy density of the Universe 
(p ~ l/(Gt 2 )) around time t ~ l/(Ga). 

Next, we consider destabilization of walls due to pressure difference Sp arising from small 
asymmetry in the conditions on the two sides. This effect competes with the two quantities 
mentioned above. Since fp ~ l/(Gt 2 ) and fx ~ (cr/(Gt 3 )) 1 / 2 , it is clear that at some point 
of time, 5p would exceed either the force due to tension or the force due to friction. For 
either of these requirements to be satisfied before to ~ V (Ga) we get 

i p >at»%.~ ( 7) 



We may read this formula by defining a factor 

where is the energy density in the wall complex immediately at the phase transition, 
which relaxes by factor T at the epoch of their decay. The factor T is strongly dependent 
on the assumed model of evolution of the wall complex, and is found to be M^/Mp t in this 
model. 



B. Evolution in matter dominated universe 



We next take up the model of evolution in which the scale factor behaves as in a matter 
dominated era by the time the domain walls get destabilised. This possibility was considered 
in [37] by Kawasaki and Takahashi. The analysis begins by assuming that the initially formed 
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wall complex in a phase tansition is expected to rapidly relax to a few walls per horizon 

volume at an epoch characterized by Hubble parameter value Hi. Thus the initial energy 

(in) 

density in the wall complex is p w ~ cHi. This epoch onwards, we assume the energy 
density of the Universe to be dominated by heavy relics or an oscillating modulus field, in 
either of which cases, the scale factor a(t) oc t 2 ^ 3 . The energy density for both of these cases 

(IT) ) 

scales as p mo d ~ P mo l I a (t) 3 '■ K the wall complex remains frustrated, i.e. its energy density 
contribution pr>w oc l/a(t), it can be seen that [37], the Hubble parameter at the epoch of 
equality of DW contribution with contribution of the rest of the matter is given by 



H eq ~ a*HfM-> , (9) 

To proceed we assume that the domain walls start decaying as soon as they dominate the 
energy density of the Universe. If the temperature at this particular epoch is To, then 
H lq ~ GT d- So from Eq.(9) we find 

Td ~ ^HfMpl (10) 
Let us assume the temperature at which the domain walls are formed T ~ a 1 / 3 . So 



4 

57T „ 4 (7 3 



From Eq.(10) we get, 



H? - Y Gas ~ (11) 



t d ~ —m ~ ~ m r hrj- ( 12 ) 



M% 2 M% 2 \M Pl 



Now requiring Sp > Tf, we get, 



3/2 



*P>mM#M (13) 



Thus in this case we find T = (Mp/Mp;) 3//2 , a milder suppression factor than in the radiation 
dominated case above.. 



C. Evolution including weak inflation 

The third possibility we consider is that the walls do not disintegrate by the time they 
come to dominate the energy density of the Universe, but in fact go on to dominate the 
energy density of the Universe. This domination however lasts for a limited epoch. Since the 
universe evolves as ait) cx t 2 , it leads to an epoch of mild inflation (as against exponential) 
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also referred to as thermal or weak inflation. This possibility has been considered [41], [42] 
in the context of removal moduli in superstring cosmology [43, 44, 45]. Such a situation is 
most likely in the case when the 5p is typically small, not large enough to destabilize the 
walls sufficiently quickly. But eventually a small pressure difference will also win over 
or fp, because either the curvature scale R diverges, as for straightened out walls, or the 
translational speed s reduces drastically. Since we have no microscopic model for deciding 
which of these is finally responsible, we introduce a temperature scale Tp at which the walls 
begin to experience instability. Note that unlike in the previous example, we will not be 
able to estimate To in terms of other mass scales and will accept it as undetermined and 
consider a few reasonable values for it for our final estimate. 

As has been studied above, at H eq the energy density of the domain wall network domi- 
nates energy density of the Universe. The scale factor at this epoch is characterized by a eq . 
Denoting the energy density of the domain walls at the time of equality as PDwiteq), the 
evolution of energy density can be written as, 

pDwitd) ~ pDwiteq) f — ) (14) 
V a d J 

where is scale factor at the epoch of decay of domain wall corresponding to time t&. If 
the domain walls decay at an epoch characterized by temperature Tp, then pDw{td) ~ Tp. 
So from Eq.(14), 

-i4 fj. \ ( a eq 



n = PDW (t eq ) i-^j (is) 

In the matter dominated era the energy density of the moduli fields scale as, 

pL d -PZ d (^-X (is) 



Substituting the value of a eq /a,d from Eq.(15) in the above equation, 

rpl2 

Ptod ~ o Ti \ ( 17 ) 

Since the energy density of the domain walls dominates the universe after the time of equality, 
PDwitd) > Pmod- So the pressure difference across the domain walls when they start decaying 
is given by, 
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where we have used the relation H^ q ~ GpDw(t eq ). Replacing the value of H eq from Eq.(9), 
and Hf ~ Gp% w ~ M|/M|,, 

^ (19) 

The JF factor introduced in eq. (8) turns out in this case to be (Tj^ Mp t ) / M}^ , rather 
sensitively dependent upon To- 



IV. SUPERGRAVITY AND LEFT-RIGHT SYMMETRY BREAKING 

The possibility that left-right symmetry may remain unbroken to low scales, and such 
breaking may be compatible with standard cosmology has been studied in an earlier work 
[29] . Specifically it was examined whether in the supersymmetric left-right symmetric model 
the parity breaking could be of hidden sector origin, and communicated to visible sector 
through gauge mediation at a low scale. The attempt is to see if several of the puzzles of the 
Standard Model and incorporation of right handed neutrinos is essentially possible within 
a few orders of magnitude of the TeV scale. It was found that messengers of a particular 
implementation of gauge mediated supersymmetry breaking can also communicate left-right 
symmetry violation. This is independent of the mechanism for the left-right symmetry 
violation in the hidden sector, the origin of which therefore remains unknown. 

The question is, does supergravity have the potential to address the origin of left-right 
symmetry violation at a low scale? There is a folk theorem that discrete symmetries can 
be violated by quantum gravity effects. The reasoning runs as follows. Formation of black 
hole horizons can cause unaccounted violation of a global charge, while preserving gauge 
charges. We then expect black hole-like virtual states in quantum gravity which can induce 
effective terms violating the global charges. Such induced terms however do not arise in 
the process of perturbative renormalization process, since every perturbative term, even in 
a non-renormalisable lagrangian should obey the expected symmetries. The effective terms 
would therefore arise from instanton-like effects. 

Due to discrete nature of the symmetry, the signal of its breaking would be in the differ- 
ence in the coefficients of the terms that get interchanged under the symmetry operation. 
The structure of supergravity ensures that at the renormalisable level gravity couples sep- 
arately to the left sector and right sector with no mixing terms. The field contents are 
identical and the gauge couplings in the two sectors are identical at this order. It appears 
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very difficult to see how supergravity would distinguish between the constants induced in 
the two sectors. We also note that N = 1 supergravity is finite at one-loop level. This leads 
us to suspect that we should not expect parity violating terms from supergravity, at least 
in the leading order in 1/Mpi. 

Thus a justification for considering 1/Mpi terms differing in their coefficients arises from 
the possibility that such are a result of gravity mediated supersymmetry breaking commu- 
nicated from the hidden sector. For this to work we must assume one of two possibilities, 
either that the gauge group governing the hidden sector does not admit left-right symmetry 
as a subgroup or that such symmetry is broken in the hidden sector. The breaking should 
then be communicated to the visible sector along with the supersymmetry breaking. The 
root cause of the parity breaking then would remain hidden as in our earlier work. We shall 
also consider the next to leading order terms and find that they may suffice only marginally 
to solve the problem of unacceptable domain walls in cosmology. 

Regardless of their origin, the structure of the symmetry breaking terms in the scalar 
potential will be similar to that of the terms that can be derived from the superpotential, 
as happens in the case of soft supersymmetry breaking terms [23]. Thus we may use the 
usual supergravity formalism to derive the terms through which Planck suppressed left-right 
symmetry breaking may get manifested. A similar approach has been adopted in the context 
of MSSM in [46] where the origin and the effect of higher dimensional operators has been 
discussed in the context of collider data. In the remainder of this section we summarize 
the essential formalism to be used in our calculation. We adopt the notation described in 
[23]. The supergravity lagrangian is obtained from three functions of complex scalar fields, 
viz. superpotential (W), Kahler potential (K) and gauge kinetic function f ab . The F-term 
contribution to the scalar potential in supergravity theory 



K 



Vf 



3e^WW*/M 2 P 



(20) 



where 



[{K- 1 )] (W* 1 + W*K l /M Pl )] 



(21) 



Making use of Eq.(21) in Eq.(20) the individual terms in Vf can be written as, 



W£WK l W*K k W l W*K k WK l 
Ml, + Ml, + Mj, 



+ Higher order terms (22) 
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Here we have considered the first term of Eq.(20). The scalar potential contains D-term 
contributions from gauge interactions which are given by 

Vd = l -Re[f~ b l D a D b ] (23) 

where 

t) a = -K\T a )U3 = —(j>* j (T a )jKi (24) 
and f ab is the gauge kinetic function which is given by 

f ab = 5 ab [l/gl + rj l /M Pl + ....] (25) 

For our purpose it will be sufficient to consider f ab = S a b/gl since we do not expect left-right 
asymmetry to arise from the gauge sector. In the following we shall consider the terms 
arising in the scalar potential from expanding W and K in the powers of (1/Mpi). 

V. GRAVITY INDUCED OPERATORS IN LEFT-RIGHT SUPERSYMMETRIC 
MODEL 

The key difference in any realistic model in comparison with the generic considerations 
of [22] is that gauge invariance restricts the structure of the non-renormalizable terms. For 
instance in [33] the operators considered were 

Kon-susv ~ ci^(A[A L ) 3 + C2 j^(A j R A R ) 3 (26) 

In other words, gauge invariance requires the terms to be 0(1/ M Pl ) rather than 0(1/ Mp{). 
A difference in pressure caused by such operators, after putting vacuum expectation values 
of the fields A ~ Mr, would be adequate to remove the domain walls accompanied by 
radiation dominated era evolution, eq. (7), and even more so the domain walls that evolved 
during an effective matter dominated era, eq. (13). 

But in the models we consider, supersymmetry forbids the kind of terms shown in eq. 
(26); on the other hand, generic parity breaking terms of 0(1/ Mpj) respecting SU(2)f and 
SU (2) r gauge invariance are permitted. This is due to additional field content with different 
gauge charges. This gain in order of 1/Mpi is however offset by the fact that due to demands 
of phenomenology, some of the fields can acquire TeV scale vacuum expectation values as 
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well. Here we have studied two successful supersymmetric implementation of left-right 
symmetry but the method can be extended to other implementations. 

The minimal left-right SUSY model is based on the gauge group SU(3) C <g> SU{2) L 
® SU(2)r (g) U(1)b-l- The anomaly free B — L global symmetry of the Standard Model 
is promoted to a gauge symmetry. The quark, lepton and Higgs fields for the minimal left- 
right SUSY model, with their respective quantum numbers under the gauge group SU(3) C 
<g> SU(2) L (g> SU(2) R ® U(1) B -l are given by, 

Q = (3,2,1,1/3), Q c = (3*,l,2,-l/3), 

L = (1,2,1,-1), L c = (1,1,2,1), 

$, = (1,2,2,0), fori = 1,2, 

A = (1,3,1,2), A = (1,3,1,-2), 

A c = (1,1,3,-2), A c = (1,1,3,2). (27) 

where we have suppressed the generation index for simplicity of notation. In the Higgs sec- 
tor, the bidoublet $ is doubled to have nonvanishing Cabbibo-Kobayashi-Maskawa matrix, 
whereas the A triplets are doubled to have anomaly cancellation. Under discrete parity 
symmetry the fields are prescribed to transform as, 

Q <- Q* c , L <-> L* c , <- $J, 

A «-> A c *, A <-> A*. (28) 

However, this minimal left-right symmetric model is unable to break parity spontaneously 
[14, 15]. Inclusion of nonrenormalizable terms gives a more realistic structure of possible 
vacua [16, 17, 47]. Such terms were studied for the case when the scale of SU(2)r breaking 
is high, close to Planck scale. 

A. The ABMRS model with a pair of triplets 

Due to the difficulties with the model discussed above, a "minimal" renormalisable model 
was developed early in [16, 17, 18], and will be referred to here as the ABMRS model. In 
this model two triplet fields Q and Q c , were added, with the following quantum numbers, 

= (1,3,1,0), ^ c = (1,1,3,0), (29) 
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which was shown to improve the situation with only the renormalisable terms [17, 18, 48]. It 
was shown that this model breaks down to minimal supersymmetric standard model (MSSM) 
at low scale. This model was studied in the context of cosmology in [28, 38] specifically, the 
mechanism for leptogenesis via Domain Walls in [49]. 

The superpotential for this model including Higgs fields only is given by, 

W LR = m A Tr AA + m A Tr A C A C + ^Tr Q 2 + ^Tr Q 2 C 
+ HijTi r 2 $f r 2 $j + aTr AQA + aTr A C Q C A C 
+ aijTr n®iT 2 <f>jT 2 + a^-Tr fl c $f ra^-ra . (30) 

Since supersymmetry is broken at a very low scale, we can employ the F and D flatness 
conditions obtained from the superpotential, to get a possible solution for the vacuum ex- 
pectation values (vev's) for the Higgs fields. 

(fi) = 0, (A) = 0, (A) = 0, 

L e 0\ (0 (A _ (0d c \ (31) 

\ o -u c ) y c oj \0 J 

This solution set is of course not unique. Since the original theory is parity invariant a 
second solution for the F and D flat conditions exists, with Left type fields' vev's exchanged 
with those of the Right type fields [38, 49]. 

With vev's as in eq. (31) the pattern of breaking is 

SU(2) L ®SU(2) R ®U(1) B „ L ^ SU(2) L ®U(1) R ®U(1) B _ L (32) 

M ^zf SU(2) L <g) U(l)y (33) 

The model introduces a new mass scale, thq. However, it was observed [17] that these terms 
can be forbidden in the superpotential by invoking an R symmetry, and then the corre- 
sponding terms appearing in the scalar potential can be interpreted as soft terms entering 
only after supersymmetry breakdown at the electroweak scale. This approach imposes a 
condition on the scales of breaking, with respect to the electroweak scale M\y, 

M R M W ~ Mi„ L (34) 
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This relation raises the interesting possibility that the scale of Mp can be as low as 10 4 
to 10 6 GeV, with corresponding very low scale 10 3 to 10 4 GeV of lepton number violation, 
opening the possibility of low energy leptogenesis [19, 49]. 

As discussed in Sec. (IV) we shall proceed to find the 1/Mpi terms in the effective potential 
by expanding Kahler potential and superpotential in powers of 1/Mpj. Here we include 
A(A), A C (A C ) and fi(fi c ) fields in the expansion of the Kahler potential and superpotential. 
The Kahler potential in this model upto (1/Mpi) can be written as, 

K = Tr (AA f + AA f ) + Tr (A c Aj + A C A\) + Tr (fiftt) + Tr (O^J) 

+ — ^(Tr AfJA 1 " + Tr AfiA t + h.c) + — -^-(Tr A c fi c Aj + Tr A c fi c A^ + /i.c). (35) 
Mpi Mpi 

Where Cl and Cp are dimensionless constants. The superpotential can also be expanded in 

powers of (1/Mpi) which is written as 



w = w ren + w n 



= m A (Tr A A + Tr A C A C ) + ^(Tr ft 2 + Tr fi 2 ) + a(Tr AfiA + Tr A c fi c A c 

+ ^_( T rAA) 2 + —^-(TrA c A c ) 2 + -^-TrA 2 Tr A 2 + —^-Tr A 2 Tr A 2 

2M Pl K } 2M Pl K cJ M Pl M Pl c 

+ ^_( T ri2 2 ) 2 + -^(Trfi 2 ) 2 + -^-M 2 TrAA + — ^— Trfi 2 TrA c A c 
4Mp; V ; 4M Pi v cJ IMpi 2Mpi 

f rri A A rn A A ^1 A 2rp A 2 k 



— TrAATrA c A c + — — TrA 2 TrA 2 + — — Tr A 2 Tr A 2 

M Pl M Pl c M Pl 

+ -^-Trfi 2 ™ 2 + -J-Trfi 2 TrA c A c + -^-Trfi 2 TrAA (36) 
Mpi Mpi Mpi 

where the coefficients appearing against the individual terms are dimensionless constants. 

The effective potential resulting from above modifications is calculated in appendix A. 

We have also examined the D-terms in the effective potential and find that they do not give 

rise to 0(^—) terms. Thus we only need to deal with F-terms. Assuming a phase in which 

the right-type fields get a non-trivial vev and all left-type fields have zero vev, the expression 

for the leading term in 1 / Mp\ the scalar potential becomes 

, rR a(cR + d R ) 4 a(a R + d R ) 3 2 
e// M~p~i R w H m~p~ M R M w \ 6 <) 

Now due to left-right symmetry, there is also a corresponding phase in which the left-type 

fields getting a vev, but no the right-type fields. For this phase the value of the effective 

potential is 

t/L a{c L + d L ) a{a L + d L ) 3 2 
e// M P ~ l R W + M~p~i MrMw (38) 
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The possibility of these two phases of approximately equal energy density gives rise to domain 
walls separating such phases. The pressure difference across the walls is proportional to the 
difference in energy density between two sides of the wall, and is given by, 

5p ~ [(c L - c R ) + (d L - d R )} f + [(a L - a R ) + (d L - d R y 



M Pl LV v " lvn M Pl 



A M R M W t A M R Myy 



,^^py (39) 

Where k a = (c^ — c R ) + (dj, — <i R ) and k m = (aj, — or) + (e^, — <i p ), and the superscript A 
refers to the ABMRS model. From Eq.(39), we see that to leading order in 1/Mpi there are 
two kinds of operators appearing in Sp, differing in powers of (Mw /Mr). 

We shall now compare these operators with the energy density required for the successful 
removal of the domain walls in the three cases to be labelled as RD, MD and WI respectively, 
discussed in Sec. (II). Comparing Eq.(7), with individual operators in Eq.(39) and taking 
the scale Mr as 10 6 GeV, and taking the more dominant term k, we get the constraint 

2 



k a > icr 10 f Mr V 

*RD>W (i 6 G evJ 



(40) 



This is easily satisfied at low scale of Mr proposed. For Mr scale tuned to 10 9 GeV needed 
to avoid gravitino problem after reheating at the end of inflation, krd ~ 10~ 4 , a reasonable 
constraint, but requires k rd to be 0(1) if the scale of Mr is an intermediate scale 10 11 GeV. 

Next, comparing Eq.(13) with individual terms in Eq.(39), the constraint on K-f /ID is found 
to be 

which puts a modest requirement on the value of krd for suitable disappearance of domain 
walls. However taking Mr ~ 10 9 GeV being the temperature scale required to have thermal 
leptogenesis without the undesirable gravitino production, leads to k md > 10 5//2 , an unac- 
ceptable requirement. The MD case is in fact generic to supersymmetric and string inspired 
models, [43, 44, 45] due to moduli production. And we find that in this case the ABMRS 
model requires a low scale of Mr and non-thermal or resonant leptogenesis. 

In the WI case, eq.(19) there is extreme sensitivity to the scales of Mr and To- Proceeding 
in same way as above comparing Eq.(19) with Eq.(39) the constraint on n is found to be, 

, 4 /10 6 GeV\ 15 / Td \ 12 

K - >10 {-Mr-) iioiev) ' ^ 

1(3 



This is a reasonable constraint for the proposed median values of the two mass scales. 
However the constraint makes the model rather strongly predictive. The scale of decay of 
the wall complex Tp can be any value below chosen Mr scale. Thus if ~ 10 4 GeV, then 
Mr is forced to be closer to the gravitino scale 10 9 Gev. This can be problematic if the 
reheating temperature after the disappearance of the domain walls in comparable to the 
temperature scale of the original phase transition. The Universe would reheat to 10 9 Gev, 
raising the possibility of unacceptable gravitino regeneration. 

Finally we consider the possibility raised in sec. IV, that the gravity induced terms are of 
direct origin, and due to one-loop fmiteness of supergravity, do not give rise to 1/Mpi terms 
in the superpotential or the Kahler potential. In such a case the most dominant operator 
to be considered is suppressed by (1/Mp;) 2 . In the ABMRS case we find, after substituting 
the vacuum expectation values, that such an operator has the magnitude 

^Next-to-leading-order ~ J\/f2 v 4c * / 

M pl 

So long as we are considering theories with Mr values less that an intermediate scale 
10 11 GeV, such terms are subdominant to the ones considered above, tightening each of 
the above constraints by a factor Mpi/Mw ~ 10 16 . Such a constraint immediately renders 
the first two scenarios of domain wall evolution cosmologically unacceptable. The third case 
of weak inflation however continues to be possible for phenomenologically acceptable values 
of the energy scales. 

To summarise the situation for the ABMRS model, we have found that there is an upper 
bound on the scale Mr if the wall evolution unfolds during a radiation dominated epoch or a 
matter dominated epoch. In the the latter case, which is generic for string theory cosmology 
with presence of heavy moduli fields, the natural value of Mr is required to be significantly 
lower than, 10 9 GeV. In the case of an evolution accompanied by weak inflationary epoch, 
there is no upper bound, rather a lower bound on the scale Mr but which is extremely 
sensitive to the value of the scale To at which the walls may finally disappear. 



B. The BM model with a single singlet 

An independent approach to improve the minimal model with introduction of a parity 
odd singlet [50], was adopted in [14, 15]. However this was shown at tree level to lead to 
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charge-breaking vacua being at a lower potential than charge-preserving vacua. 

Recently, an alternative to this has been considered in [51] here a superfield S = (1, 1, 1, 0) 
also singlet under parity is included in addition to the minimal set of Higgs required as 
in Eq.(27). The A c , A c are required for SU(2) R <g> U(1)b-l symmetry breaking without 
inducing R-parity violating couplings. The singlet field S is introduced so that SU(2)r 
(g) U(1)b~l symmetry breaking occurs in the supersymmetric limit. The superpotential 
including the Higgs fields is given by, 

The superpotential is given by, 

W LR = + 

Where 

W w = h; i} L T r 2 <l> t r 2 L c + hf Q T r 2 <S> t T 2 Q c 
+ ii*L T T 2 AL + HL cT t 2 A c L c 

+ S [ X* Tr AA + A Tr A C A C + \' ab Tr ^r 2 ^ b r 2 - M 2 R ] (44) 
W {2) = M A Tr AA + M* A Tr A C A C 

+ /i ab Tr ^T 2 $ b T 2 + M S S 2 + A S ,S 3 (45) 

In this analysis the terms in has been assumed to be zero. Dropping the terms in 
makes the theory more symmetric and more predictive. It is observed that dropping 
quadratic and cubic terms in S leads to an enhanced i?-symmetry. Further, dropping the 
massive couplings introduced for A's means that A masses arise purely from SUSY breaking 
effects, keeping these fields light and relevant to collider phenomenology. Dropping the 
fi a b terms for $ fields makes it possible to explain the \x parameter of MSSM as being 
spontaneously induced from S vev through terms in W^ 1 ). Additionally, absence of the 
terms can be shown to solve the SUSY CP and strong CP problems. 

The presence of linear terms in S in makes possible the following SUSY vacuum, 

(S) = 0, Xv R v R + X*v L v L = M 2 R (46) 

where v L (v L ) and v R (v R ) are the vev's of the neutral components of A(A) and A C (A C ) fields 
respectively. From Eq.(46) it is clear that we have a flat direction in the vl - v R space. 
Assuming that the flat directions are lifted, we have two choices viz. 

vr = v r = 0, \v L \ = \v L \ = -j= (47) 

V A* 
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vl = v L = 0, |wr| = = — ^= (48) 

The important result is that after SUSY breaking and emergence of SUSY breaking 
soft terms, integrating out heavy sleptons modifies the vacuum structure due to Coleman- 
Weinberg type one-loop terms which must be treated to be of same order as the other terms 
in V e ff . Accordingly, it is shown[51] that the V e ^ contains terms of the form 

VflL^c) ~ -lifmicTr (A c At)Af - |/| 2 mi c Tr (A c A c )Tr (AjAj)^ (49) 

where Af and A2 are constants obtained from expansion of the effective potential. Pres- 
ence of these terms is shown to lead to the consequence of a preference for electric charge- 
preserving vacuum over the charge-breaking vacuum, provided m| c < 0. 

Further Eq.(47) also constitutes a valid solution of Eq.(46). In this vacuum the soft terms 
can give rise to the following terms in the effective potential, 

*£&»(A) ~ "I/I Vl* (AAt)Af - |/| 2 m 2 Tr (AA)Tr (A* A*) A*; (50) 

with A\ and A\ constants. Thus the choice of known phenomenology is only one of two 
possible local choices, and formation of domain walls is inevitable. 

Here we analyse the full superpotential without setting terms zero. Like in the 

ABMRS model, we study here the scalar potential taking supergravity into account. In- 
cluding triplet fields and singlet field S, the Kahler potential in this model upto 0{\/M P i) 
is given by, 

K = Tr (AA 1 " + AA 1 ") + Tr (A C A* + A C A\) + \S\ 2 + -^(Tr ASA 1 " + Tr AS At + h.c) 

Mpi 

+^(Tr A C SA\ + Tr A C SA\ + h.c) + -^(S 3 + h.c). (51) 

Mpi Mpi 

The superpotential upto (1/Mpi) order is given by, 
W = W ren + W nr 

= m A (Tr AA + Tr A C A C ) + M S S 2 + X S S 3 + S[X*Ti AA + ATr A C A C - Af|] 

+ — (TrAA) 2 + — ^— (TrA c A c ) 2 + -^-TrA 2 Tr A 2 + — TrA 2 Tr A 2 

2Mpi ' 2Mpi ' IMpi 2Mpi 

+ + —S 2 Tr AA + -^S 2 TrA c A c 

M Pl M Pl M Pl 

+ J-Tz AATr A C A C + -^Tr A 2 Tr A 2 + -^Tr A 2 Tr A 2 (52) 
Mpi Mpi Mpi 
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The effective potential has been calculated in the appendix(B) considering the term from the 
effective potential (K~ 1 )* k l Wj^W l . When the right-type fields get a vev the scalar potential 
can be written as, 

v 5f ~ ^~Ml + Jg-sM* + ^LsV R + (53) 
Mpi Mpi Mpi Mpi 

where s is the scale at which S gets a vev. To calculate the potential when only the left 
type fields get vacuum expectation values, we introduce corresponding coefficients a^, etc. 
We then compute the pressure difference across the walls as 

^ M R t„ „ n S * M R 

5p ~ (a L - a R )—^- + . . . + (C L - C R ) 



M Pl y " llJ M Pl 

.,bs 3 M 2 r 



Ml lB s 3 M 2 R 
-^ + ... + k' B —-^ (54) 



M Pl M Pl 

where k b = {a>L — ap), k' b = {Cl — Cp), superscript B referring to the BM model and 
... is in lieu of terms which, as we explain next, are relatively unimportant. According to 
BM model, the value s is of the scale of supersymmetry breaking. If this scale is TeV scale, 
then Mp is expected to be higher and the n' B terms are subdominant. In case however the 
supersymmetry breaking scale is 10 11 GeV, it could be higher then the scale of Mp. In this 
case the k' b term is expected to dominate. 

In the case of TeV scale supersymmetry breaking, comparing the first term in Eq.(54) 
with Eqs.(7), (13) and (19) of with sec. (II), we get the corresponding constraints on possible 
values of k as 

^ > (lo^v) • < 55 > 



» > in- (^s\ 



1/2 



"MD 



Vio 6 Gevy 

M R 



Thus for the proposed Mp scale of 10 6 GeV there is no serious constraint on k b values. 
Only in the scenario with weak inflation, if Tp> scale is high, such as lOOGeV, the value 
Mp ~ 10 6 GeV becomes marginal, but due to large powers of the mass scale present, a small 
increase in Mp easily offsets the effect of the increase in Tp>. Overall, the kind of operators 
obtained in this particular model provide no constraint on the mass scale Mp, as long as 
the scale of supersymmetry breaking is TeV scale. 
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To check other possibilities, we consider the supersymmetry breaking scale Ms and hence 
s to be ~ 10 11 GeV. In this case we get, proceeding as above, 

12 /-,^fi^_Tr\ 17 /mllnTA 3 



/ Mg \ ^ /KFGeVV 

Vio 6 Gevy V s y 

K - >10 llOGevJ l^rj l^^J (60) 



This shows that there is no particular constraint on the induced parity breaking coefficients 
due to increase in the scale of supersymmetry breaking Ms, so long as Mr < s = Ms ~ 
10 n GeV. 

In summary the BM model remains mostly unrestricted by the present considerations. 
This is due to newer terms possible with a gauge singlet. 



VI. CONCLUSION 



Ever since the discovery of massive neutrinos it has become a tantalising possibility that 
the small neutrino masses arise from rich physics at a high energy scale, which in turn 
would naturally incorporate right handed neutrinos. A model that treats this new content 
symmetrically with the known contents would naturally lead to the requirement of parity 
symmetry in the high energy model. If this discrete symmetry is spontaneously broken it 
would lead to formation of domain walls in the early Universe. 

We have considered three scenarios for the evolution of transient domain wall networks 
ending in their decay. We characterize each model by a dimensionless parameter T which 
is the ratio of the available pressure difference across a wall at the time of its decay to the 
the characteristic energy density M^ in the Universe at the time of formation of the wall 
complex, see Eqs. (7), (13), (19). Of the three scenarios, the first one unfolds entirely 
in a radiation dominated universe, in which the dynamics is governed by the interplay of 
forces due to friction and tension and the pressure difference across the walls. Here we find 
that the parameter T is given by (Mr/Mpi) 2 . The second scenario unfolds in a matter 
dominated era, where the domain wall complex decays as soon as the energy density of the 
same dominates the energy density of the Universe. The parameter T in this case is given 
by {M R /Mpif 2 . The third scenario is an extension of the second one where we assume 
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that the domain wall network comes to dominate the energy density of the Universe and 
continues to do so for a finite epoch before it decays. Characterizing the required pressure 
difference Sp across the domain walls in this case requires additional input, the ratio of 
scale factor values a eq /ad where subscript eq refers to the epoch at which the domain walls 
become equally as important as the rest of the matter and subscript d refers to the epoch 
at which the decay of the domain walls occurs. We characterize this ratio by an equivalent 
"decay temperature" To defined in Eq. (15). The formula for T in this case shows very 
high sensitivity to the mass scales concerned. 

For each of these cases we study the viability of the two specific models of spontaneous 
parity breaking, the ABMRS model, sec. V A and the BM model, VB. Each of the particle 
physics models permits intrinsic operators that must match up to the required parameter J 7 , 
resulting in final destablisation of the wall complex. If the operators available in the model 
cannot provide a 5p of required magnitude, the wall complex would not be destabilized, 
leading to unacceptable cosmology. 

The ABMRS model turns out to be more restrictive, using as it does only non-trivial 
representations of the gauge group. In this case high scale of parity breaking becomes 
conditionally disfavored, though still viable if the wall evolution leads to a weak inflationary 
epoch. The BM model containing a singlet turns to not be restricted by the considerations 
here. 

Our main conclusion is that a low scale scenario with Mr ~ 10 6 GeV or lower is viable 
and generic. Specifically, in the ABMRS model with domain wall evolution in a matter 
dominated epoch Mr is restricted to remain less than 10 8 GeV, eq. (41). A matter dominated 
epoch is generic to string theory inspired models with occurance of moduli fields of mass 
scale 10 9 GeV and hence this restriction is of special interest. 

As a broader conclusion we learn that gravity induced higher dimensional operators 
can effective in ensuring the removal domain walls generic to a theory with spontaneous 
breakdown of parity, but the result can be model dependent. This mechanism is then an 
alternative to an earlier studied [29] possibility of parity breaking mediated by the messengers 
in a version of gauge mediated supersymmetry breaking scenario. The constraints on the 
parameters for that scenario to ensure disappearance of domain walls were rather stringent. 
We now see that supergravity induced terms can ensure wall disappearance with a very 
modest constraints, which can be predictive in some cases. The question of what underlying 
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physics results in gravity inducing parity breaking terms remains open. It appears also tied 
to whether the terms arise directly from supergravity in the visible sector or through the 
gravity mediated supersymmetry breaking effects. 
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APPENDIX A: AMS MODEL 



The scalar potential contains D-term contributions from gauge interactions. From Eq.(23) 

V D = ^Re(g 2 J5 ab D a D b ) 
= ^Re(D a D a ) 

= ^-Re\Tr(K A (T a )A) + Tr(K A (T a )A) + Tr(K n (T a )tt)\ 2 
a 2 

= ^-Re\Tr(K A r a A) + Tr(K A r a A) + Tr{K n r a n)\ 2 
8 

+ ^-Re\Tr(K A A) - Tr(K A A)\ 2 
a 

= ^-Re\2Tr(A^r a A) + — -Tr(ftAV a A) 
8 Mpi 

2C 

+2Tr(AV a A) + — -Tr(fiAVA) 

Mpi 

2C 

+4Tr(Or a fi) + — -[Tr(A^Ar a n) + Tr(A^ Ar a n)]\ 2 
Mpi 

a 12 

+^- J Re|2Tr(A t A) - 2Tr(A t A) 
8 



Tr(AfiAt) - -±Tr{A^)\ 2 (Al) 

Mpi Mpi 

The D-term vanishes after putting the vev's for the corresponding fields. From above it is 
clear that we can not find 1/Mpi suppressed terms from Vb- So we have to go for Vp to find 
the desired terms. Here we consider the first term appearing in Eq.(22) i.e. (k" 1 )* W^W . 

Substituting the Eqs.(36) and (35) in Eq.(22), the terms which contribute are: 

V R ~ ^Trn 2 Jin c A c Al + ^-m A TrA c A c TrA c At 
+ ^m A TrA c A c TrA c At + ^M c 2 M c A c At 

IVlpi AM pi 

+ — mJM^ + — -Trfi^AcficAc + terms higher order in l/M P i (A2) 
Mpi Mpi 
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In the ABMRS model we have the relation: 

M|_ L ~ M R M W ; u = -1^1 =M R ;d = d= (2m A m n /a 2 ) 1/2 = M B „ L (A3) 

After putting the vev's for the corresponding fields and making use of appropriate scale, the 
terms upto 0(1/Mpi) are 

Vr ^ ajo^^ + aja^^ 

Mpi Mpi 
APPENDIX B: BM MODEL 

The leading order terms in this model comes from the first term of Eq.(22). writing 
explicitly the individual terms 

v * = ( K ^)a c aI W l wAc + i*- 1 ) KM W L wAc + (K^ss* W sW s (Bl) 

Calculating the above terms the terms in the scalar potential in lowest order in 1/Mpi are 
given by, 

V F rsj -^m A TrA c A c TrA c AJ + -^A*,STr A c A c TrA c AJ 
Mpi c M Pl 

+ [|A| 2 5 3 TrA c A]; + M A STrA c Aj + S 2 (A*M A TrA c Aj + AM A TrA c Aj)l 
Mpi 

+ -^-AS , TrA c A c TrA c AJ + ^M s ^ 2 TrA c A c + -^s 3 TrA c A c 
Mpi c M Pl M Pl 

+ Other terms (B2) 

After putting the vev's for the neutral components of the triplet field and using the appro- 
priate scale the term in the highest power of Mr, 

Vr ~ — Mn + ^sM R + —s 2 M 3 R + —s 3 M 2 + Other terms (B3) 
M Pl M Pl R M Pl Mpi K ' 
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